function [ phi ] = read_phi( n_kt,beta )

    phi = zeros(size(n_kt,1),size(n_kt,2));
    for k = 1:size(n_kt,1)
        for t = 1:size(n_kt,2)
            top = sum(n_kt(k,t))+beta(t);
            db = 0;
            for v = 1:size(n_kt,2)
                db = db + n_kt(k,v)+beta(v);
            end
            phi(k,t) = top/db;
        end
    end

end

